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ABSTRACT 

Detection of 21 cm emission of H I from the epoch of reionization, at redshifts z > 6, is limited primar¬ 
ily by foreground emission. We investigate the signatures of wide-held measurements and an all-sky 
foreground model using the delay spectrum technique that maps the measurements to foreground 
object locations through signal delays between antenna pairs. We demonstrate interferometric mea¬ 
surements are inherently sensitive to all scales, including the largest angular scales, owing to the nature 
of wide-held measurements. These wide-held effects are generic to all observations but antenna shapes 
impact their amplitudes substantially. A dish-shaped antenna yields the most desirable features from 
a foreground contamination viewpoint, relative to a dipole or a phased array. Comparing data from 
recent Murchison Wideheld Array observations, we demonstrate that the foreground signatures that 
have the largest impact on the H I signal arise from power received far away from the primary held of 
view. We identify diffuse emission near the horizon as a signihcant contributing factor, even on wide 
antenna spacings that usually represent structures on small scales. For signals entering through the 
primary held of view, compact emission dominates the foreground contamination. These two mecha¬ 
nisms imprint a characteristic pitchfork signature on the “foreground wedge” in Fourier delay space. 

Based on these results, we propose that selective down-weighting of data based on antenna spacing 
and time can mitigate foreground contamination substantially by a factor ^100 with negligible loss 
of sensitivity. 

Keywords: cosmology: observations — dark ages, reionization, hrst stars — large-scale structure of 
universe — methods: statistical — radio continuum: galaxies — techniques: interfero¬ 
metric 
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1. INTRODUCTION 

At the end of the recombination epoch, the Universe 
was completely neutral. This period, referred to as the 
Dark Ages in the Universe’s history, is characterized by 
the localized accumulation of matter under the influence 
of gravity. And it ended with the formation of the first 
stars and galaxies which started emitting ultra-violet and 
X-ray radiation, thereby reionizing the neutral medium 
in their surroundings. This commenced the epoch of 
reionization (EoR) - a period of nonlinear growth of mat¬ 
ter density perturbations and astrophysical evolution. 


Observing redshifted 21 cm radiation generated by 
the spin flip transition of H I has been identifi ed as 


a direct probe of the EoR (Sunvaev & Zeldovich 

1972; 

Scott & Rees 19901: iMadau et al. 19971: iTozzi et al. 

2000; 

Iliev et al.l I2002D. Detecting this signal has recentlv 


emerged as a very promising experiment to fill the gaps 
in our understanding of the Universe’s history. 


Sensitive instruments such as the Square Kilometre Ar¬ 
ray (SKA) are required for direct observation and to¬ 
mography of redshifted H I. Numerous pathfinders and 
precursors to the S KA such as the Mu r chison Wide- 
held Array (MWA: iLonsdale et al.l 120091: iTingav et al.l 
120131: lBowman~~et al. 20 131), the Low Frequency Ar¬ 
ray (LOFAR; Ivan Haarlem et all l2013h . and the Pre¬ 
cision Arr ay for Pr obi ng th e Epoch of Reionization 
(PAPER: IParsons et al.l f20Toh have become operational 
with enough sensitivity for a st atistical detect ion of 
the EoR H I power spectrum (Bowma n et al.l [ 2 006; 
Parsons et al. 2012a; Beardsley et al. 2013; Dillon et al.l 
20131: iThvagaraian et al.l 12013c iPober et al. 203). The 
Hydrogen Epoch of Reionization Arra\P^I (HERA) is cur¬ 
rently under construction using new insights gained with 
the MWA and PAPER. 

A key challenge in the statistical detection of the red¬ 
shifted H I 21 cm signal, via the spatial power spec¬ 
trum of temperature fluctuations, arises from the con¬ 
taminatio n by Galactic and extragalactic foregrounds 
(see, e.g. . iDi Matteo et al. 2002; Zaldarria ga et al.ll2004 
Furlanetto et al . 120061: lAli et al.l 120081: fe ernardi et al.l 

2009. 120101: IGhosh et aljf2012E iMorales Hewittl (12004) 

show that the inherent isotropy and symmetry of the 
EoR signal in frequency and spatial wavenumber ( k ) 
space make it distinguishable from sources of contami¬ 
nation which are isolated to certain k mo des b y virt ue of 
their inherent spectral smoothness ([Morales et al.ll 2 006; 

Bowman et a l. 2009; Liu & Tegmark 2 0 lit IParsons et ahl 

2012bl : iDillon et al.l 120131: IPober et al.ll2013l ). Since this 
contamination is expected to be several orders of magni¬ 
tude stronger than the underlying EoR H I signal, it is 
critical to characterize foregrounds precisely in order to 
reduce their impact on EoR H I power spectrum detec¬ 
tion sensitivity. 


Considerable effort is being made toward understand¬ 
ing the k- space behavior of foreground signatures in the 
observed power spectrum and form ulating robus t esti- 


mators of the true power spectrum (Bowman et al. 

2009; 

Liu et al. 20091 Datta et al. 2010; Liu & Tegmark 

2011; 

Morales et al.li2012fc ITrott et al. 20121 

3 ober et al. 

2013; 

Thvagaraian et al.l 120131: IDillon et al. 

2014 |Liu et al.l 

2014alb|). A model that provides a generic explana- 


22 http://reionization.org/ 


tion for the observed foreground power spectrum has 
emerged, whereby the wide-held (and chromatic) re¬ 
sponse of the instrument causes the power in smooth 
spectrum foregrounds to occupy higher k- modes into the 
so-called “wedge”. The conservative foreground strategy, 
referred to as avoidance , that has developed alongside 
this work is to discar d k- modes which could be contam¬ 
inated (e.g., IParsons et alJ : 2014 ). The more aggressive 
alternative is to subtract a sky model and regain access 
to modes that would be discarded by avoidance. In both 
cases, which parts of the sky are most critical to either 
avoid or subtract has remained largely uncertain. Here, 
we focus primarily on extending the avoidance strategy 
by identifying foreground components at greatest risk to 
“leak” from foreground modes to EoR modes and propos¬ 
ing a scheme for down-weighting these components. 

Foregrounds with intrinsic deviations from spectral 
smoothness, instruments with high chromaticity, polar¬ 
ization leakage, calibration errors, or approximations in 
power spectrum analyses can contaminate the true EoR 
H I power spectrum. Here we use existing catalogs 
and a high fidelity instrumental model to capture both 
foreground and instrumental chromaticity. To decouple 
these effects fr om possi b le analysis e ffects, such as those 
pointed out bv lHazelton et al l ion . we comp ute powe r 
spectra using a per-baseline approach of IParsons et al.l 
(2012b). This approximates the power spectrum as the 
inverse Fourier transform of the spectra generated by the 
instrument’s correlator. 

In $2] we provide an overview of the delay spectrum 
technique. We investigate signatures generic to all wide- 
held measurements of EoR power spectra in $3] In $H 
we present the foreground model and a variety of instru¬ 
ment models to rank antenna shapes based on foreground 
contamination. In $5j we describe the MWA setup, sum¬ 
marize the observing parameters, and present the result¬ 
ing data. Simulations using these observing parameters 
are compared with the data and analyzed for foreground 
signatures. We report two important findings: fore¬ 
grounds that most severely obscure the redshifted 21 cm 
power spectrum are not caused by emission in the cen¬ 
tral held of view, but rather by bright objects from near 
the horizon; and, diffuse Galactic emission plays a sig¬ 
nificant role hitherto unpredicted. In $6l we offer an ini¬ 
tial description of a more precise foreground avoidance 
technique that minimizes foreground contamination us¬ 
ing prior knowledge of the sky to down-weight adversely 
contaminated baselines. We present a summary of our 
work and hndings in $71 

2. DELAY SPECTRUM 

We prov ide a short overview of the delay spectrum 
technique (jParsons et al.l 12012 alibi ). 

Interferometer array data known as visibilities , U&(/), 
represent correlations between time-series of electric 
helds measured by different antenna pairs with separa¬ 
tion vectors b and then Fourier transformed along the 
time axis to obtain a spectrum along the frequency (/) 
axis. If /(s,/) and A(s,/) are the sky brightness and 
antenna’s directional power pattern, respectively, at dif¬ 
ferent frequencies as a function of direction on the sky 
denoted by the unit vector (s), and W[(f) denotes in¬ 
strumental bandpass weights, then V&(/) can be written 
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as: 

v b (f) = jj A(a, f) I(a, f) Wi(f) e~ i2 ^ dfi, ( 1 ) 

sky 


it is closely related to the sought power spectrum con¬ 
taining critical information about spatial scales. In using 
the delay spectrum technique, visibilities from different 
baselines have not been averaged together. 


where, c is the speed of light, and dQ is the solid an¬ 
gle element to which s is the unit normal vector. This 
equation is valid in general, including wide-held mea¬ 
surem e nts, and is a s light adaptation from Ivan Cittertl 
(|1934fl . IZernikel ([1938D , and [Thompson et al.l (12001D . 

The delay spectrum , Vb(r ), is defined as the inverse 
Fourier transform of Vb(f) along the frequency coordi¬ 
nate: 

V b (r) = J V b (f ) W(f) e i2 ^ T d/, ( 2 ) 

where, W (/) is a spectral weighting function which can 
be chosen to control the quality of the delay spe ctrum 
(|Vedantham et al.l 120121: IThvagaraian et al.l 120131) . and 
r represents the signal delay between antenna pairs: 


The delay transform conventions used in this paper are 
described in appendix [A] Vb(r) is expressed in observer’s 
units of Jy Hz. 

The delay spectrum has a close resemblance to cosmo¬ 
logical H I spatial power spectrum. Appendix [B] gives an 
overview of the similarities and differences between the 
two. Foregrounds can be described in either framework. 
For our study, we find the delay spectrum approach to 
be simple and yet extremely useful. 

In order to express a quantity derived from Vb(r) whose 
units are the same as that of the cosmological H I power 
spectrum, we define the delay power spectrum: 


P d (fcL,fc||) = \V b (r)\ 2 


(DaA 

\X 2 AB J 


(d 2 ad\ 

\~AB~) 


with 


k± a 


2A!) 

D 


_ 2nr f 2 iH 0 E(z) 
fc|l= c(l + *)2 



( 4 ) 

( 5 ) 

(6) 


where, A e is the effective area of the antenna, A B is the 
bandwidth, A is the wavelength of the band center, k b is 
the Boltzmann constant, /21 is the rest frame frequency 
of the 21 cm spin flip transition of H I, z is the redshift, 
D = D(z ) is the transverse comoving distance, AD is 
the comoving depth along the line of sight corresponding 
to AH, and h, Hq and E{z) = [£2m(1 + ^) 3 + ^k(l + z) 2 + 
^ A ] 1/2 are standard terms in cosmology. Throughout the 
paper, we use (2 m = 0.27, = 0.73, £2k = 1 — CIm ~ ^a, 

Hq = 100 km s -1 Mpc -1 , and -Pd (&_!_> &||) is in units of 
K 2 (Mpc/h) 3 . 

In summary, the delay spectrum, V 5 (r), is obtained 
from visibilities which are the basic data blocks measured 
by each antenna pair, using equations Q] and [ 2 ] 14 (t) cap¬ 
tures all the effects of EoR H I signal corruption caused 
by foregrounds and the instrument. At the same time, 


2 . 1 . Delay Space 

We give a brief overview of some parameters of Fourier 
space which are generic to all experiments that use a sim¬ 
ilar approach. Figure [l] illustrates the Fourier space in 
which the delay (and power) spectra of redshifted H I 
observations are calculated. | 6 | and fcj_, denoting spa¬ 
tial scales in the transverse direction (tangent plane to 
the celestial sphere), form the x-axis. r and /cy, de¬ 
noting spatial scales along line of sight form the y- axis. 
Foreground emission maps to a wedge-shaped region in 
Fourier space, hereafter referred to as the foreground 
wedge dDatta et al.l [ 2010 lh whose boundaries are deter¬ 
mined by the antenna spacings and the light travel times 
across them. These boundaries, called horizon d elay 
limits ([Vedantham et al1l2 01 2t [Parsons et al.ll2012bt h are 
shown by solid lines. 


k ± [h Mpc 1 ] 

0.031 0.063 0.094 0.126 



Figure 1. Fourier space in which delay (and power) spectra of 
EoR H I signals are calculated. The x-axis is denoted by \b\ (an¬ 
tenna spacing) or k± (transverse wavenumber). The y-axis denoted 
by t (delay) or foy (line of sight wavenumber). Here, kj_ and k y are 
obtained for a frequency of 185 MHz. The dark shaded region is 
referred to as the foreground wedge where smooth spectrum fore¬ 
grounds reside. Its boundaries (solid lines), given by light travel 
time for corresponding antenna spacings, are referred to as hori¬ 
zon delay limits. Narrow extensions of the wedge (white unshaded 
strips) are caused by convolution with the instrument’s spectral 
transfer function. Regions excluding the wedge are expected to 
be relatively free of foreground contamination and are generally 
referred to as the EoR window. There are undesirable grating 
responses (dotted-dashed lines) specific to the MWA. Hence, we 
conservatively identify a restricted region of high EoR sensitivity 
(medium shade) and refer to it as the MWA EoR window. 

The spectral transfer function of the instrument con¬ 
volves the foreground wedge and stretches it further (un¬ 
shaded narro w strips bounded by so li d and dashed lines) 
along r-axis (jParsons et al.l l2012bl : IThvagaraian et all 
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I2013D . The width of this narrow strip is inversely pro¬ 
portional to the operating bandwidth. The region of 
Fourier space excluding the foreground wedge and the 
narrow strips is the so-called EoR window , shown in light 
and medium shades. In the context of EoR studies in 
Fourier space, the H I power spectrum from the EoR is 
expected to decrease rapidly with \k\. Hence, the bright¬ 
est EoR signal will be observed on the shortest baselines 
and smallest delays. Thus the regions of interest for EoR 
studies relying on avoidance strategy are just beyond the 
horizon delay limits (dashed lines) on short baselines, 
marked as regions of maximal EoR sensitivity. 

In the specific case of the MWA, which has a passband 
constructed using coarse channels, there are period grat¬ 
ing responses resulting in repetitions of the foreground 
wedge at multiples of 0.78 /as. Thus the MWA EoR win¬ 
dow lies outside the dashed lines but inside the first grat¬ 
ing response (|r| < 0.78/is, dotted-dashed lines) and is 
shown in medium shade. 


2.2. Delay Spectrum Deconvolution 

We obtain the delay spectrum of visibilities by taking 
the delay transform of each baseline’s spectrum (equa¬ 
tion choosing W(f) to be a Blackman-Harris win¬ 
dow function. The sky spectrum is multiplied in the 
instrument by the instrumental passband and flagging 
of frequency channels possibly corrupted by radio fre¬ 
quency interference, which together are represented by 
the weights, W[(f). In delay-space, these weights trans¬ 
late into a convolution by a point spread function (PSF). 
We deconvol ve thi s PSF using a one dimensional CLEAN 
algorith m (iTavlor et all 1 199911 as described for the de- 
lay axis (jParsons fc Backed 120091: iParsons et al.l l2012bl ) 
to obtain the final delay spectra. The CLEAN procedure 
iteratively finds and subtracts peak values convolved by 
the Fourier transform of the weights. We limit the se¬ 
lection of peaks to modes inside the horizon delay limit, 
corresponding to smooth spectrum objects in the visible 
sky hemisphere. 


3. WIDE-FIELD MEASUREMENTS 



sky 


where, Z, m, and n denote the direction cosines to¬ 
ward east, north, and zenith respectively, with n = 
y/1 — l 2 — m 2 , and: 


d^ 


dZ dm 

\/l — l 2 — rn 2 


(8) 


When the synthesized field is small, where A(s, f ) or 
7(3,/) is significant only for |Z| 1 and \m\ 1, equa¬ 

tion [7 reducestoashmleJwoMimen^iomlFomier trans¬ 
form ( Tavlor et al.lll999t [Thompson et ahll2Q01[) between 
the apparent sky brightness and measured visibilities. It 
is in this context that radio interferometers are under¬ 
stood to be sensitive only to fluctuations and not to a 
uniform sky brightness distribution. 

In a wide-held measurement, neither A(s, /) nor 
7(3,/), in general, is negligible anywhere in the visible 


hemisphere. The solid angle per pixel on the sky in direc¬ 
tion cosine coordinates changes significantly with direc¬ 
tion (equation [8]) , increasing steeply toward the horizon. 
Hence, the approximations in the narrow-held scenario 
do not apply. For example, even if A(s,f) and 7(3,/) 
are held constant across the visible hemisphere, the am¬ 
plitude of the integrand in equation [7] is still dependent 
on direction. Therefore, in a signihcant departure from a 
narrow-held measurement, the wide-held visibility from 
a uniform brightness distribution on a non-zero antenna 
spacing is not zero. 

Figure [2] shows the wide-held delay power spectrum of 
a uniformly illuminated sky with no spectral variation 
as measured by antenna elements arranged identical to 
that in the MWA antenna array layout (jBeardslev et all 
l2012f ) with a uniform power pattern across the sky and a 
bandwidth of 30.72 MHz centered around 185 MHz (re¬ 
fer to 34.11 for a detailed description of the instrument 
model). Notice the steep rise in power toward the hori¬ 
zon limits. These wide-held effects are prevalent on all 
antenna spacings, including the longest ones used in this 
study. 


k ± [h Mpc 1 ] 


0.031 0.063 0.094 0.126 k 2 (Mpc/h) 3 



50 100 150 200 


|b| [m] 

Figure 2. Wide-field effects on delay power spectra produced 
with a uniform sky brightness distribution measured by antenna 
pairs with a uniform power pattern across the visible hemisphere. 
Delay power spectra are obtained using equations □ n and 0] for 
each baseline, which are then stacked by baseline length. The axes 
correspond to cosmological dimensions. The non-zero response of 
the interferometer array to a uniform brightness distribution and 
the prominent edge brightening close to the horizon delay limits are 
wide-held effects. These are prevalent on all antenna spacings and 
are generic to all instruments used in wide-held measurements. 


We interpret this as due to equal-sized delay bins sub¬ 
tending larger solid angles near the horizon thereby con¬ 
taining larger integrated emission. Further, baseline vec¬ 
tors (including those with largest lengths) are foreshort¬ 
ened toward the horizon along their orientation. Thus, 
they become sensitive to larger angular scales that match 
the inverse of their foreshortened lengths along these di¬ 
rections. 
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iThvagaraian et ah (120131) found evidence of this fea¬ 
ture in their statistical models. In line with our rea¬ 
soning, they attribute it to a steep rise in solid an¬ 
gles su btended by delay bins near the horizon limits. 
iPober et al.l (|2013f ) also find a similar “edge brightening” 
feature which they attribute to Galactic plane emission 
near the horizon. From their discussion, it is unclear 
what fraction of power in that feature arises from such 
wide-held effects. 

We conclude these are generic to all instruments mak¬ 
ing wide-held measurements. The nature of the specihc 
instrument used for observing will control the amplitude 
of these effects, which we explore below. 

4. SIMULATIONS 

We describe the instrument and foreground models 
used in our simulations. 

4.1. Instrument Model 

In our present study, we use a latitude of -26?701 
and an antenna layou t identical to that of the MWA 
(Beardsl ey et al.l f2012h for the observatory. The array 
is arranged as a centrally condensed core of ^ 300 m — 
there are many spacings in the range 5-50 m — and a ra¬ 
dial density that falls oh as the inverse of the radius, with 
the longest baselines at 3 km. Here we focus on antenna 
spacings \b\ < 200 m (spatial scales relevant to reion¬ 
ization). Their deviation from coplanarity is negligible. 
For geometrical intuition, we restrict the orientation (# 5 , 
measured anti-clockwise from east) of all baselines to lie 
in the range —67?5 < 0b < 112?5. Baselines oriented in 
the other half-plane measure conjugate visibilities with 
delays of equal magnitude but of opposite sign and hence 
are ignored in our analysis. We choose an observing fre¬ 
quency of 185 MHz ( z ~ 6 . 68 ) and a hat passband of 
width A B = 30.72 MHz to roughly match those of ongo¬ 
ing MWA EoR observations (discussed in detail in $5]). 

One of the principal components of the instrument 
model is the antenna power pattern, A(s, /) (see equa¬ 
tion!]]). It is determined by the shape of its aperture. Us¬ 
ing a few examples, we examine the role the geometrical 
shape of the aperture plays in shaping the characteris¬ 
tics of delay power spectrum. We consider the following 
antenna elements placed at the MWA tile locations: 

1. Dipole: an east-west dipole of length 0.74 m at a 
height 0.3 m above a ground plane. A e = (A/ 2 ) 2 . 

2 . Phased Array : a 4x4 array of isotropic radiators 
with a grid spacing of 1.1 m at a height 0.3 m above 
the ground plane placed in an arrangement similar 
to that of an MWA tile. A e = 16 (A/ 2 ) 2 . 

3. Dish: diameter of 14 m similar to that proposed 
for HERA, with A e « 154 m 2 . The power pat¬ 
tern is simulated using an Airy pattern where its 
sensitivity beyond the horizon is forced to zero. 

4.2. Foreground Model 

In wide-field measurements, it is important to consider 
an all-sky model for foreground objects in evaluating the 
features seen in the power spectrum instead of restrict¬ 
ing only to the primary field of view, a point also sup¬ 
ported by Pober et al. (2015, in preparation). We use a 


foreground model that includes both diffuse and bright 
compact components. 

For the diffus e com p onent, we use an all-sk y radio fore¬ 
ground model (jde Oliveira-Costa et al.ll2008l ) to estimate 
the emission at 185 MHz. At this frequency, since this 
map is predominantly based on the 408 MHz map of 
lHaslam et al.l (|1982l) which has an angular resolution of 
0?85, we smoothed the 185 MHz map to the same res¬ 
olution. However, to avoid any artifacts from sampling 
this map, we sample it at ~ 27' intervals. We model the 
diffuse foreground spectra with a unique spectral index 
at each pixel in the map, estimated from model maps at 
170 MHz and 200 MHz. 

The model described above is primarily a model 
of the diffuse foreground sky. While it contains 
faint compact emission blended in with the diffuse 
e missio n, brigh t po int sou rces have been removed 
(jde Oliveira-Costa et al.l [20081 ). In order to supplement 
it with missing bright compact emission, we use clas¬ 
sical radio source confusion estimates to determine the 
nominal flux density threshold and include point sources 
brighter than this threshold. Slightly different criteria 
are in common use in radio astronomy to estimate ra¬ 
dio so urce confusion (see Appendix of lThvagaraian et al.l 
I2013L and references therein). For an angular resolution 
of 0?85, using a conservative ‘ S c = 5cr c ’ criterion, we de¬ 
termine the flux density threshold to be ~ 10 Jy. Other 
liberal criteria that yield a lower threshold carry a greater 
risk of double-counting point sources which might be al¬ 
ready blended in with the diffuse sky model. 

We u se a combination of the NRAO VLA Sky Survey 
(NVSS; iCondon et al.|[l998l) at 1.4 GHz and the Sydney 
University Molonglo Sk y Survey (SUMSS; iBock et all 
119991: iMauch et al.ll2003[) at 843 MHz to provide our point 
source catalog due to their complementary survey foot¬ 
prints covering the entire sky, and matched flux density 
sensitivity and angular resolution. The SUMSS catalog 
covers the sky with declination S < —30° with a lim¬ 
iting peak brightness of 6-10 mJy/beam and an angu¬ 
lar resolution of ^ 45". The NVSS covers the sky with 
S > —40° with a similar angular resolution and a limit¬ 
ing flux density of ~ 2.5 mJy for point sources. 

From the SUMSS catalog, we select objects whose de¬ 
convolved major axes are equal to 0", thereby strictly 
selecting point sources. From the NVSS catalog, we ex¬ 
cluded objects that overlap with those in the SUMSS 
survey footprint. Point sources from NVSS were selected 
if the convolved major axes were not greater than « 47", 
which matches the angular resolution of the survey. Us¬ 
ing a mean spectral index of (q s n ) = —0.83 (flux density, 
S(f) oc f asp ) obtained by iMauch et al.l (2003) for both 
NVSS and SUMSS catalog objects, we calculate the cor¬ 
responding flux densities at 185 MHz, S'iss- From this 
subset, we choose point sources with S 185 > 10 Jy. The 
selection of such bright point sources is not affected by 
minor differences in flux density sensitivity of the two 
surveys. We verified that our selection criteria ensure a 
similar areal density of objects in the two surveys. 

These criteria yield 100 objects from the SUMSS cat¬ 
alog and 250 objects from the NVSS catalog. Together 
with the diffuse foreground model, we obtain an all-sky 
foreground model consisting of both compact and diffuse 
emission. Figure [3] shows the diffuse (top) and compact 
(bottom) foreground emission model used in our study. 
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In this snapshot pointed toward zenith at 0.09 hr LST, 
the Galactic center in the diffuse model has just set in 
the west. 
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Figure 3. Foreground model at 185 MHz consisting of diffuse 
emission (top) in units of K and bright point sources (bottom) in 
units of Jy, visible during a snapshot at 0.09 hr LST. In the diffuse 
model, the Galactic center has just set in the west. Color scales 
are logarithmic. 


4.3. Role of Antenna Geometry 

The power patterns of the aforementioned antenna ge¬ 
ometries at 185 MHz for this zenith pointing are shown 
in Figure Iffal 

The delay power spectra without thermal noise com¬ 
ponent for these antenna shapes are shown in Figure l4bl 
The occupancy of the power patterns on the sky is clearly 
correlated with that in the delay spectra. Further, the 
strength of the primary lobe centered on the pointing 
center is correlated with the delay power spectrum cen¬ 
tered on r = 0; and, the overall rate of decrease in the 
power sensitivity away from the pointing center is corre¬ 
lated with the rate of drop in power away from r = 0. 

The levels of foreground contamination in the EoR 
window varies substantially across the different antenna 
shapes: ~ 10 4 K 2 (Mpc/h) 3 , < 10 2 K 2 (Mpc/h) 3 , and 
< 1 K 2 (Mpc/h) 3 for the dipole, phased array, and dish, 
respectively. The severity of foreground contamination 
inside the foreground wedge both in strength and occu¬ 
pancy also evidently decreases as the antenna element 
is changed from a dipole to a phased array to a dish. 
For instance, notice that the foreground contamination 
in k- modes between fcy =0 and the horizon limits de¬ 
creases from 10 5 K 2 (Mpc /h) 3 in a phased array to 
^ 10 K 2 (Mpc /h) 3 in a dish. As a consequence, k- modes 
in the foreground wedge that may be deemed too con¬ 
taminated for EoR studies in the case of a dipole or a 
phased array can potentially become accessible when us¬ 
ing a dish. 


Finally, a distinct feature common to all these aper¬ 
ture shapes is that the foreground contamination near 
the horizon delay edges is significant even on wide an¬ 
tenna spacings (> 10 5 K 2 (Mpc/h) 3 ). We have argued 
this arises due to wide-field effects. The prevalence of 
this feature across different antenna shapes demonstrates 
it is generic to all wide-field measurements. The ampli¬ 
tude of this effect, however, can be controlled via choice 
of antenna shape and through weighting of aperture illu¬ 
mination. A dish-shaped antenna appears to hold a sig¬ 
nificant advantage over a dipole or a phased array from 
the viewpoint of foreground contamination. 

Typically, the sensitivity of antennas to the primary 
field of view is high compared to the rest of regions on 
the sky. Combined with the wide-field effects seen ear¬ 
lier, it leads to a “pitchfork”-shaped signature inside the 
foreground wedge , as exemplified in the case of a dish. Al¬ 
though the exact appearance of this signature depends 
on the antenna power pattern, we use the term pitchfork 
hereafter, to broadly refer to the combination of fore¬ 
ground power in the primary field of view and the en¬ 
hancement of foreground power near the horizon limits 
due to the nature of wide-field measurements. 

5. THE MURCHISON WIDEFIELD ARRAY 

We now use our simulations to analyze features in ob¬ 
served dela y power spectrum ob t ained using the M WA 
instrument ([Lonsdale et al.ll2009b iTingav et al1l2013l) . 

MWA construction was completed in 2012 and, af¬ 
ter commissioning, began its EoR observing program in 
2013. The MWA is a 128-tile interferometer capable of 
observing a 30.72 MHz instantaneous band anywhere in 
the range 80-300 MHz. Each tile is a phased array of 
16 dipoles, each in the shape of a bow-tie. This yields a 
primary field of view >20° wide and multiple secondary 
lobes. See Bear dsley et all d2012) for the tile layout. 

The MWA passband of width A B = 30.72 MHz is di¬ 
vided coarsely into 24x1.28 MHz sub-bands with each 
sub-band weighted by a digital filter. The coarse chan¬ 
nel shape is obtained using an eight-tap polyphase filter 
bank (PFB) and a Kaiser window with parameter /3 m 5. 
Each of these coarse bands consists of 32 fine channels 
of width 40 kHz each. After correcting for the shape of 
these coarse channels, the fine channels at the edges of 
these sub-bands are flagged because they are known to 
be contaminated by aliasing at a low level. 

The MWA is expected to be sensitive to the power 
spectrum of the H I signal over the redshift range 
6 < z < 10 ([B owma n et all 120061 : iBeardslev et al.l 120131: 
iThvagaraian et al.ll2013[ ). Over 600 hr have been cur¬ 
rently observed usi ng the MW A, targe ting science objec¬ 
tives outlined in lBowman et al. (120131 ). 

The MWA targets two primary low-foreground fields 
for reionization observations. Here, we focus on the field 
at R.A. = 0 h , deck = —30°. The MWA tracks a patch 
of sky through antenna beams formed and steered elec¬ 
tronically by controlling delay settings of the dipoles in 
an MWA tile. The pointing system is optimized to points 
on a regular ^7° grid. The sky drifts across to the near¬ 
est available pointing, shifting between grid points (once 
every ^ 30 minutes). This process is repeated through¬ 
out the course of the observation « 4.86 hr. 

The observations used here were acquired on 2013 Au¬ 
gust 23. We have chosen two sections of duration 112 s 
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(a) Power patterns 
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(b) Delay power spectra 

Figure 4. Power patterns (top panels) and simulated delay power spectra (bottom panels) for different antenn a sh apes at 185 MHz 
centered on zenith. Antenna shapes used are: dipole (left), phased array (middle), and dish (right). Refer to H4.ll for details of the 
antenna models. The strength and occupancy of the power patterns are correlated with those of delay power spectra. White dotted lines 
in the delay power spectra mark the boundaries of the foreground wedge determined by the horizon delay limit and antenna spacing. The 
foreground wedge and the EoR window are most severely contaminated in the case of the dipole while it is the least for the dish. The 
phased array has intermediate levels of contamination. Foreground power close to the horizon delay limits in all three cases is significant 
even on long baselines. The foreshortening of baselines toward the horizon makes them sensitive to foreground emission on large size scales. 
The amplitude of this feature strongly depends on the shape of the antenna element. It is highest for a dipole (which has a strong response 
near the horizon) and least for a dish. 


each from this night for detailed study. These were cho¬ 
sen to provide a selection of possible foreground and in¬ 
strumental conditions. As an example of a nominal ob¬ 
serving setup we choose a zenith pointing; as an example 
of poor foreground conditions, we choose a pointing when 
the field is ~ 2 hr from zenith. This pointing has a signif¬ 
icantly higher secondary lobe structure and is observed 
when the bright galactic center is well above the horizon. 

These two pointings are at LST 22.08 hr and 0.09 hr, 
which are hereafter denoted as off-zenith and zenith 
pointings, respectively. 

5.1. Initial Data Processing 

Th e data are flagged for interference (jOffringa et al.l 
EqU), removing 3% of the data and averaged in time and 
frequency from the raw 0.5 s, 40 kHz to 2s, 80 kHz. These 
data are then calibrated to a simulation of the sky con¬ 
taining 2420 point-like objects selected from the MWA 
Commissioning Survey (MWACS; iHurlev-Walker et~ahl 
1201 4 ). It has a flux density limit of 25 mJy and a declina¬ 


tion range of —12° to —40° evenly covering the field of 
view of the observations reported here. The objects used 
in calibration are selected to lie inside the 5% contour of 
the primary lobe of the tile power pattern. The calibra¬ 
tion algorithm — base d on forward modeling software 
bv I Sullivan et al.l (2012) an d the calibration method de¬ 
scribed bv lSalvini fe Wiinholdsl (|2014t ) — computes com¬ 
plex gain solutions per channel per antenna averaged to 
two minute intervals. The solutions are fairly low signal 
to noise so we iteratively average along the antenna and 
frequency dimensions to capture the relatively indepen¬ 
dent passband and antenna-to-antenna variation. First, 
we average the channel gains over all antennas to obtain a 
high signal to noise measurement of the bandpass. After 
applying this single passband, we do a second round of 
calibration and fit second and first order polynomials for 
amplitude and phase respectively for each antenna. This 
flattens any residual variation in bandpass and removes 
small phase slopes due to variations in cable delay. Fi¬ 
nally, we fit for an additional phase known to be caused 
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by small reflections in a subset of cables. 


5.2. Modeling 

The MWA tile power pattern is modeled as a mutually- 
coupled 4-by-4 dipole array with the overall power pat¬ 
tern of each individual dipole cal culated via finite el¬ 
ement electromagnetic simulation (jSutinio et al.ll20l4f ). 
To speed up simulations, we find that a phased array of 
isotropic radiators at a height of 0.3 m above an infinite 
ground plane provides a very good approximation to the 
full simulation, hence we use the idealized dipoles. We 
also assume that each individual dipole signal has ran¬ 
dom delay fluctuations of rms 0.05 ns, a number in line 
with the known rep eatability and stabil ity level of the 
analog signal chain ([Bowman et al.ll2007l) . Besides hav¬ 
ing the effect of adding a time-dependent uncertainty in 
the power pattern, these random delay fluctuations re¬ 
duce the coherence in the phased addition of dipole sig¬ 
nals resulting in deviations from predicted models of the 
power pattern, most prominently at its nulls. 

We use the model described in 94.21 for the foreground 
sky. Figure [5] shows the diffuse emission and bright point 
source foreground models for the two chosen pointings 
with the modeled MWA tile power pattern contours over¬ 
laid. Notice the presence of a portion of the Galactic 
plane and the bright Galactic center in the westward sky 
in the diffuse sky model, where the MWA tile power gain 
is significant (> 12%). In the zenith pointing, the Galac¬ 
tic plane has set and the power pattern in that direction 
is at least 16 times smaller. 

We estimate thermal noise, AV , in the observed visi¬ 
bilities, Vb(f), using the rms of V&(r) obtained from data 
after delay-deconvolution across all antenna spacings for 
\r\ > 1 /is using the relations: 


AV = ^N^AVAf, 


AV 


2 T sys 

AeV2 A/At’ 


and 


(9) 

( 10 ) 


where, A/ = 80 kHz, At = 112 s, and N c h = AH/A/ is 
the number of frequency channels. The choice of thresh¬ 
old for r is well outside the foreground window , where 
foreground contamination is negligible and thus yields a 
robust estimate of ^~sys- We find the average system tem¬ 
perature to be ^ 95 K. Hence, for our simulations, we 
use ^~sys = 95 K to match the thermal noise in data. 


5.3. Comparison of Data and Model 

With the aforementioned foreground model, and in¬ 
strumental and observational parameters, we simulate 
visibilities using equation [lj Figure [6] shows the delay 
power spectra from off-zenith and zenith pointings ob¬ 
tained from MWA observations and modeling. Notice 
the qualitative agreement of amplitude and structure be¬ 
tween the two. The Galactic center and the Galactic 
plane visible in the off-zenith pointing make it appear 
brighter in the foreground wedge as a branch with r < 0. 

In order to make a quantitative comparison of delay 
spectra obtained with MWA data and our simulations, 
we consider the uncertainty in the assumed spectral in¬ 
dex of our foreground model. Our foreground models 
are derived from other higher frequency catalogs and 
sky maps. The inherent spread in spectral index in¬ 


creases the uncertainty while predicting fluxes at the ob¬ 
serving frequency. Using simple error propagation, the 
fractional error in the delay spectrum amplitude caused 
by the spread in spectral index is ln(/orig//) Aa SP 5 
where, / or i g is the original frequency at which the cat¬ 
alog or map was created, / = 185 MHz is the MWA 
observing frequency, and Aa sn is t he spread (HWHM) 
in spectral index. From iMauch et al. (120031) . we assume 
Aa sp « 0.35 for point sources f rom NYSS and SUMSS 
catalo gs. Although the model of Ide Oliveira-Costa et all 
(120081 yields a spectral index per direction on the sky, 
we could assume similar uncertainties exist in spectral 
indices of our diffuse sky model as well, wh ich is predo m¬ 
i nantl y derived from the 408 MHz map of lHaslam et al.l 
([1982D . Thus, fractional errors in delay spectrum ampli¬ 
tudes from compact and diffuse components are ^70% 
and ~30% respectively. 

In addition to intrinsic model uncertainty, delay spec¬ 
tra from simulations and data each have fluctuations 
due to thermal noise in the delay spectrum with rms 
~ 1.4xl0 7 Jy Hz. We estimate the ratio of delay spec¬ 
tra from data and simulations as p = |U 6 d (t)| / | U 6 S (r) |, 
where superscripts D and S denote data and simulation, 
respectively. The median absolute deviation of log 10 p 
inside the foreground wedge for both pointings is « 0.28. 
This corresponds to ~ 90% fractional difference between 
data and modeling on average with either pointing. 

We also simulated delay spectra after assigning spec¬ 
tral indices drawn randomly from a gaussian distribu¬ 
tion with a mean of (<a sp ) = —0.83 and a HWHM of 
Att sp = 0.35 to the point sources in our compact fore¬ 
ground model. These simulations typically yielded a me¬ 
dian absolute deviation of ~ 0.29 for log 10 p indicating 
fractional differences of ~ 95% between different realiza¬ 
tions. This demonstrates that a fractional deviation of 
~ 90% observed between data and simulations is in line 
with expectations when the aforementioned uncertainty 
in foreground models, thermal noise fluctuations in mea¬ 
surements, and uncertainties in antenna power pattern 
due to random delay fluctuations are taken into account. 

These uncertainties are presented only to confirm the 
qualitative agreement already seen between data and 
modeling in Figure [6] These estimates are conserva¬ 
tive. A full treatment of all uncertainties and devia¬ 
tions from ideal behavior such as fre quency dependent er¬ 
rors i n tile power patte rn ([Bernardi et al. 2015), calibra¬ 
tion (Da tta et al.ll2010l) . data corruption due to interfer¬ 
ence, anisop lan atic wide-held imaging and ionospheric ef¬ 
fects (jlntema et all 200 9) will bring the simulations much 
closer in agreement with observations, but is beyond the 
scope of this paper. Hereafter, our focus is to explore in 
detail the foreground signatures embedded in the fore¬ 
ground wedge of the MWA instrument. 

5.4. Analysis of Foreground Signatures 

Having shown that the simulation matches the data to 
the level of expected uncertainties, we proceed to exam¬ 
ine in further detail the key signatures seen in simulated 
delay spectra. A number of factors are responsible for 
the characteristics noted in the delay spectra obtained 
from data and through simulations. In subsequent sec¬ 
tions, we provide a detailed explanation of our results 
as a combination of these factors. Note that numerous 
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(a) Diffuse Sky Model (b) Bright Point Source Model 

Figure 5. Sky brightness temperature of the diffuse foreground model (left) and flux densities of bright point sources (right) at 185 MHz 
visible during off-zenith (top) and zenith (bottom) pointings. The color scales are logarithmic. MWA tile power pattern contours are 
overlaid. The contour levels shown are 0.00195, 0.00781, 0.0312, 0.125, and 0.5. The Galactic center and a portion of the Galactic plane 
are prominently visible during the off-zenith pointing in the diffuse sky model and the MWA tile power gain is significant (> 12%) in that 
direction. In contrast, emission from the Galactic plane in zenith pointing is significantly less. 


features may overlap at different degrees of significance 
depending on combinations of parameters. We assign 
the features to their predominant causes. Secondly, we 
have used noiseless cases to clearly illustrate the observed 
foreground signatures. With the addition of noise in the 
visibilities, some of the weaker features may not be as 
prominently visible. Since the foreground signatures are 
far too numerous and subject to a multitude of parame¬ 
ters like baseline length and orientation, power pattern, 
patch of sky under observation, and instrumental config¬ 
uration, we highlight only the most notable features in 
the foreground delay power spectra. 

Figure [7| shows the delay power spectra obtained from 
the diffuse (left) and compact (right) foreground emission 
for the off-zenith (top) and zenith (bottom) pointings 
without thermal noise component. Some of the notable 
signatures are discussed below. 

5.4.1. Galactic Center on Eastward Antenna Spacings 

The most prominent signature seen in the off-zenith 
pointing (top left panel, Figure [7j) is due to the bright 
Galactic center situated on the western horizon co¬ 
located with one of the bright secondary lobes of the 
power pattern. It appears as a bright branch near 
the negative delay horizon delay limit. This feature is 
strongest at short antenna spacings and fades with in¬ 
creasing antenna spacing. The bright signature is absent 
in the zenith pointing (bottom left panel, Figure [7j) be¬ 
cause the Galactic center has set below the horizon. 

5.4.2. Ubiquitous Diffuse Emission 

Diffuse emission outside the Galactic plane manifests 
in the primary field of view as a branch at r > 0 and 
r = 0 in the off-zenith and zenith pointings respectively. 
The former is seen at r > 0 because the primary lobe of 


the power pattern is centered eastward of zenith, whereas 
in the latter it is centered at zenith. As we see from Equa¬ 
tion U each baseline measures a single spatial mode on 
the sky with an angular size scale inversely proportional 
to the length of the baseline projected in the direction 
of the emission. Thus, in the zenith pointing, the hor¬ 
izontal line at r = 0, fades away on antenna spacings 
| 6 | > 125 m because the diffuse sky model is devoid of 
spatial structures on scales < 0?75. 

5.4.3. Diffuse Emission on Wide Antenna Spacings 

In both pointings the diffuse emission (left panels) is 
prominent near the horizon delay limits extending to the 
widest antenna spacings. This is a characteristic signa¬ 
ture of the wide-held effects discussed in JU It is evident 
at all LSTs in our simulations. Thus, diffuse emission 
from far off-axis directions manifests as an edge-heavy 
two-pronged fork across all baselines. It decreases in 
strength with increasing baseline length but is neverthe¬ 
less present in all baseline orientations. 

5.4.4. Compact Foreground Signatures 

In contrast to the delay power spectra of diffuse emis¬ 
sion, compact emission (right panels) manifests as a 
center-heavy structure in either pointing. 

The amplitude response of an interferometer to a point 
source is, to first order, hat across baseline length. Since 
the primary held of view in the off-zenith pointing is cen¬ 
tered eastward of zenith, the bulk of the compact fore¬ 
ground emission is seen in a branch with r > 0. In the 
zenith pointing, compact emission from the same patch 
of sky is seen as a bright horizontal arm at r = 0 since 
the primary held of view is centered at zenith. 

Foreground emission at r = 0 and r < 0 in the off- 
zenith pointing is caused by point sources co-located with 
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Figure 6. Delay power spectra from MWA data (left) and modeling (right) for the off-zenith (top) and zenith (bottom) pointings. The 
foreground wedge is bounded by white dotted lines . Model matches the data to a level consistent with the uncertainties in foreground 
models and the antenna beam (as discussed in 


secondary lobes of the power pattern. On the other hand, 
point sources co-located with secondary lobes of power 
pattern in the zenith pointing are revealed as faint but 
distinct branches at positive and negative delays depend¬ 
ing on the orientation of antenna spacing and direction 
of emission on the sky. 

5.5. The “Pitchfork” 

Delay spectra from the foreground model in our study 
display a composite feature set drawn from the features 
of compact and diffuse foreground models. Here we com¬ 
pare the relative strengths of emission from different spa¬ 
tial scales in our composite foreground model. 

When not dominated by the bright emission from the 
Galactic center, the delay power spectrum of the com¬ 
bined foreground model is composed of diffuse and com¬ 
pact emission both of which are significant. This is illus¬ 
trated by a detailed examination of the zenith pointing. 

Figure [8] shows delay power spectra of three antenna 
pairs of different spacings oriented northward during the 
zenith pointing; each is a different vertical slice of the 
delay power spectra plots shown in Figure [3 The dif¬ 
fuse, compact, and composite components are shown as 
solid red, cyan, and black lines, respectively. The horizon 
delay limits are shown as a pair of vertical dotted lines. 

The peak at r = 0 (corresponding to the primary 
lobe in the power pattern) with a value of ^ 10 7 - 
10 8 K 2 (Mpc/h) 3 , independent of antenna spacing, is pre¬ 
dominantly determined by compact emission. The peak 
at r = 0 from diffuse emission is ^ 10 3 times fainter and 
decreases rapidly with increase in antenna spacing. This 
is the response expected from different antenna spacings 
toward compact and diffuse emission. 

Near the horizon delay limits, the diffuse component is 
brighter relative to the compact component. Here, dif¬ 
fuse emission does not decrease as rapidly with increas¬ 


ing antenna spacing as was seen at r = 0. In fact, even 
on widely spaced antennas, diffuse emission in the delay 
power spectrum near the horizon delay limits exceeds 
that in the primary lobe by about t hree o rders of magni¬ 
tude. This feature is described in £15.4.31 and attributed 
to wide-held measurement effects discussed in JU 

Simulations with the complete foreground model show 
the combination of center-heavy features dominated by 
compact emission in primary held of view, and edge- 
heavy features from both types of emission especially 
the diffuse component near the horizon. This results 
in a characteristic pitchfork structure imprinted in the 
foreground wedge and should be evident in observations. 

The observability of the pitchfork signature predicted 
in this paper depends on the relative levels of uncertainty 
in the foreground model and fluctuations from thermal 
noise. In our simulations, since thermal noise in these 
very short duration snapshots is rsj 10 4 K 2 (Mpc/h) 3 and 
features near the horizon delay limits are also of compa¬ 
rable amplitudes, the pitchfork feature is not expected 
to be detected, though this feature is marginally visible 
in the zenith pointing of observed data (see Figure [6]). 
We attribute this to differences between our foreground 
model and the actual sky. Deeper observations should 
reveal the feature clearly. 

We also note that increasing the antenna spacing pro¬ 
gressively improves the resolution along the delay axis by 
increasing the number of delay bins inside the foreground 
wedge. This improves the localization of foreground ob¬ 
jects whose signatures are imprinted in the delay power 
spectrum. For instance, there is an increase in the num¬ 
ber of secondary peaks in the delay power spectrum be¬ 
tween r = 0 and horizon delay limits as the antenna 
spacing increases from ~84 m to ~1T1 m. In this case, 
these correspond to secondary lobes of the power pattern 
that he between the primary lobe and the horizon along 
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Figure 7. Simulated delay power spectra (in units of K 2 (Mpc/h) 3 ) for for the diffuse (left) and compact (right) foreground models in the 
off-zenith (top) and zenith (bottom) pointings without any thermal noise. The axes and color scale are identical to those in Figure [6] In 
the off-zenith pointing, emission from the Galactic center is the most prominent feature seen as a branch at r < 0. In the zenith pointing, 
delay power spectrum from diffuse emission has a two-pronged fork -shaped structure and is present even at wide antenna spacings due to 
wide-held effects. Compact emission is centrally concentrated. 


the local meridian. At short antenna spacings, due to 
relatively lower resolution along the delay axis inside the 
foreground wedge and a consequent loss of localization 
of foreground emission, these secondary peaks blend in 
with other major peaks and are not distinctly visible. 

6. baseline-based foreground mitigation 

Here, we investigate the susceptibility of particular an¬ 
tenna spacings to foreground contamination arising out 
of bright foreground objects located near the horizon and 
present a technique to substantially mitigate such con¬ 
tamination. We use the MWA as an example. 

The Galactic center in the off-zenith pointing is one 
such example already available in our study. Figure l9al 
shows the sky model (top: compact component, bottom: 
diffuse component) in this pointing. The Galactic center 
is the most dominant source of foreground contamina¬ 
tion from the diffuse sky model and is co-located with 
a bright secondary lobe of the power pattern near the 
western horizon. Figure [9b] shows the sky mapped to de¬ 
lays registered by the baseline vectors, of length 100 m 
for instance, oriented toward north (top panel) and east 
(bottom panel). Figure [9cl shows the delay spectra on 
baselines oriented northward (67?5 < #b < 112?5) at 
the top and eastward (—22?5 < #b < 22?5) at the bot¬ 
tom. The Galactic center manifests itself most distinctly 
near the negative horizon delay limit on short eastward 
baselines in the delay power spectrum (bottom panel of 
Figure [9c]). Consequently, the spillover caused by the in¬ 
strument’s spectral transfer function from the foreground 
wedge into the EoR window affects the northward base¬ 
lines the least and is most severe on eastward baselines 
(particularly the short ones) evident by the bright verti¬ 
cal stripes of foreground contamination. 

With a foreground model known a priori in which 


structures and locations of very bright foreground ob¬ 
jects such as the Galactic center or AGN are available, 
we can predict the response across antenna spacings as a 
function of observing parameters such as LST, power pat¬ 
tern, etc. This allows us to programmatically screen data 
for antenna spacings that are severely contaminated by 
foregrounds near the horizon delay limits. These can be 
weighted appropriately during data analysis. We demon¬ 
strate such a screening technique, whereby we use the 
bright object’s location and structure to discard antenna 
spacings of certain lengths and orientations to mitigate 
foreground contamination in the EoR window. 

In our example, we discard eastward antenna spacings 
(—22?5 < #b < 22?5) of lengths |b| < 30 m. Foreground 
contamination was found to be insensitive to removal of 
wider antenna spacings as discussed below. Figure fTTTI 
shows the delay spectra obtained with all antenna spac¬ 
ings (top panel) and after applying our screening tech¬ 
nique (bottom panel) on the off-zenith observation. No¬ 
tice the significant reduction in foreground spillover into 
the EoR window via the removal of bright vertical stripes 
on short eastward antenna spacings. 

This screening technique can be generalized to opti¬ 
mize between foreground mitigation and loss of sensi¬ 
tivity from discarding data. Figure [TT] shows how the 
typical foreground contamination^ in the MWA EoR 
window depends on the orientations and lengths of dis¬ 
carded antenna spacings. We choose antenna spacings 
oriented eastward to varying degrees of directedness, i.e., 
—7?5 < #b < 7?5 (solid circles), —15° < #b < 15° (solid 
squares), and —22?5 < #b < 22?5 (solid stars). Among 
antenna spacings that satisfy these criteria, we discard 

23 Foreground contamination is measured by standard deviation 
of noiseless Pd(k±, k n) from foregrounds in the MWA EoR window. 
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Figure 8. Simulated delay power spectra for three chosen north¬ 
ward oriented antenna spacings of length: ~84 m (top), ~104 m 
(middle), and ~171 m (bottom). The baseline length and orienta¬ 
tion is specified in each panel. The solid red, cyan, and black lines 
denote contributions from diffuse, compact, and composite fore¬ 
ground models respectively. Vertical dotted lines mark the horizon 
delay limits. Compact emission dominates the central regions of 
the delay power spectra while both components, especially diffuse 
emission on short antenna spacings, dominate near the horizon de¬ 
lay limits, giving rise to a characteristic pitchfork -shaped structure. 

data from those whose lengths are shorter than |6| max 
(x-axis) and show foreground contamination estimated 
in the EoR window from all remaining antenna spacings. 

In other words, Figure [Tl] demonstrates the progress 
in foreground mitigation as orientation and maximum 
length of discarded antenna spacings are varied. The 
fraction of discarded antenna spacings discarded relative 
to the total number is shown in dotted lines for different 
ranges of 6% It is seen that foreground contamination 
can be mitigated by a factor between ^ 2 (|#b| < 7?5) 
and ~ 100 (|#b| < 22?5). The latter limit is achieved 
with a mere 5% loss of data for |6| max — 30 m. Dis¬ 
carding antenna spacings with lengths | 6 | > 30 m does 
not mitigate foreground contamination any further and 
would only lead to loss of sensitivity as the fraction of 
discarded baselines increases from ~ 5% to ^ 25%. 

If there was a bright point source in place of the Galac¬ 
tic center, it will give rise to foreground contamination 
even on longer antenna spacings. Such cases will ne¬ 
cessitate discarding more or all of the eastward antenna 
spacings. The MWA was used as an example. In gen¬ 
eral, the direction, strength and type of emission, and 
the array layout will determine such thresholds. 


In principle, instead of discarding selected antenna 
spacings altogether, we could down-weight them based 
on an optimal scheme. For instance, the estimates of 
covariance computed from the delay transform bins can 
be naturally fed into the cov ariance-we ighted power spec¬ 
trum estimation techniques (jLiu et al .1112014allbf ) . It could 
also be used to downweight or flag contaminated base¬ 
lines in imaging applications. This technique provides a 
very simple and yet effective tool to suppress the effects 
of foreground contamination in EoR data analysis. 

7. SUMMARY 

Our primary motivation in this work is to understand 
how the various bright foregrounds will manifest in three- 
dimensional power spectrum of H I from 21 cm reion¬ 
ization observations. In units of temperature variance, 
the dynamic range between bright foregrounds and the 
21 cm signal is expected to be ^ 10 s ; a detailed un¬ 
derstanding of how foregrounds can corrupt the 21 cm 
power spectrum is therefore essential. This analysis ex¬ 
tends previous work by simulating the entire sky rather 
than just the central field of view and by providing a 
comparison with early observations with the MWA. By 
making use of the delay spectrum technique to estimate 
the power spectrum, we are able to observe the effects 
of foregrounds while avoiding entanglements with more 
sophisticated power spectrum estimators. 

We find that all wide-held instruments, typical of 
modern EoR observatories, imprint a characteristic two¬ 
pronged fork signature in delay spectra. This arises from 
two related effects: delay bins near the horizon subtend 
larger solid angles and therefore contain larger integrated 
emission; and, foreshortening of baselines toward the 
horizon makes them sensitive to emission on large angu¬ 
lar scales which match the inverse of their foreshortened 
lengths. These effects combined with higher sensitivity 
of antennas in the primary held of view results in a char¬ 
acteristic pitchfork signature. The amplitude of these 
generic signatures can be controlled by careful design of 
antenna aperture. In contrast to a dipole and a phased 
array such as an MWA tile, a dish such as the one pro¬ 
posed for HERA is found to yield the least foreground 
contamination and thus preferable for EoR studies. 

Simulating in many important respects the response of 
the MWA to an all-sky foreground model that consists 
of diffuse Galactic emission and bright point sources, we 
conhrm that the modeled delay spectra are in agreement 
with data obtained with the MWA to within expected 
uncertainties in foreground models. 

Our simulations enable us to identify numerous signa¬ 
tures of different components of foreground emission seen 
in the delay spectra. We establish the relationship be¬ 
tween these signatures and observing parameters such as 
antenna pointing and LST, instrument parameters such 
as antenna power pattern, and foreground parameters 
such as type of emission, etc. 

The bright Galactic center at the edge of the western 
horizon co-located with one of the far secondary lobes 
of MWA tile power pattern is the brightest source of 
foreground contamination in the off-zenith pointing. It 
manifests itself near the negative horizon delay limit in 
the delay power spectrum on eastward antenna spacings. 

Diffuse emission in the primary field of view is promi¬ 
nent on shorter antenna spacings. However, it is also 
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Figure 9. (a) Sky model showing compact (top) and diffuse ( bottom ) emission (adopted from Figure [5]). The Galactic center is very 

prominent in diffuse emission on the west co-located with a bright secondary lobe of the power pattern, (b) Sky hemisphere mapped to 
delays observed on antenna spacings oriented north (top) and east (bottom). Delays vary linearly with antenna spacing length. Color scale 
shown is for a 100 m antenna spacing. The bright Galactic center will appear at r = 0 in northward antenna spacings and close to negative 
horizon delay limit on eastward antenna spacings. (c) Simulated delay power spectra on antenna spacings oriented northward (top) and 
eastward (bottom). White lines denote horizon delay limits. The bright Galactic center is prominently visible close to negative horizon 
delay limit, especially on short eastward antenna spacings. These are also the most severely contaminated by foreground spillover. The 
northward antenna spacings, on the other hand, are the least contaminated. 
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Figure 10. Simulated delay spectra power for the off-zenith point¬ 
ing with all antenna spacings included (top) and with short east¬ 
ward antenna spacings discarded (bottom). Discarded antenna 
spacings (black vertical stripes) have lengths \b\ < 30 m (leftward 
of vertical dashed line) and orientations |#b| < 22?5. The spillover 
from the bright Galactic center near the negative horizon delay 
limit from the foreground wedge is lowered by a factor ~ 100 when 
short eastward antenna spacings are discarded. 


prominent near the horizon limits even on wide antenna 
spacings — an effect of the wide-held nature of the 
measurement. On the other hand, compact emission 
predominantly maps onto central regions of the fore¬ 
ground wedge. Features arising from compact emission 
co-located with primary and secondary lobes of the an¬ 
tenna power pattern have been identified. In general, de- 


Figure 11. Drop in foreground contamination in the MWA EoR 
window , and loss of data for the off-zenith pointing as a function 
of discarded baselines. Eastward baselines with varying degrees of 
directedness— |#b| < 7? 5 (solid circles), |6>b| < 15° (solid squares), 
and |#b| < 22? 5 (solid stars) — and lengths |b| < |fo| max (:c-axis) 
are discarded. Loss of data (dotted lines) is measured by discarded 
baselines as a fraction of the total number for the corresponding 
cases. Foreground contamination in the EoR window (solid lines) 
drops by a factor ~2 (|6>b| < 7? 5) to ~100 (|6>b| < 22? 5). The latter 
limit can be achieved with a mere 5% loss of data at |6| max — 30 m, 
and discarding longer baselines (\b\ > 30 m) has no effect in further 
reducing foreground contamination. 

lay power spectrum signatures of compact emission are 
center-heavy while those of diffuse emission are edge- 
heavy which results in a characteristic pitchfork signa¬ 
ture. This will be evident when thermal noise is suffi- 
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ciently lowered, as larger volumes of data are processed. 

We also provide a simple and effective tool based on 
the delay spectrum technique that can potentially mit¬ 
igate foreground contamination by nearly two orders of 
magnitude in EoR data analysis by discarding or down¬ 
weighting data from antenna pairs most affected by fore¬ 
ground contamination, with negligible loss of sensitivity. 
In conclusion, we find that inclusion of emission models, 
both diffuse and compact, all the way to the horizon is 
essential to explaining the observed power spectrum. 


This work was supported by the U. S. National Science 
Foundation (NSF) through award AST-1109257. DCJ is 
supported by an NSF Astronomy and Astrophysics Post¬ 
doctoral Fellowship under award AST-1401708. JCP is 
supported by an NSF Astronomy and Astrophysics Fel¬ 
lowship under award AST-1302774. This work makes 
use of the Murchison Radio-astronomy Observatory, op¬ 
erated by CSIRO. We acknowledge the Wajarri Yamatji 
people as the traditional owners of the Observatory site. 
Support for the MWA comes from the NSF (awards: 
AST-0457585, PHY-0835713, CAREER-0847753, and 
AST-0908884), the Australian Research Council (LIEF 
grants LE0775621 and LE0882938), the U.S. Air Force 
Office of Scientific Research (grant FA9550-0510247), 
and the Centre for All-sky Astrophysics (an Australian 
Research Council Centre of Excellence funded by grant 
CE 110001020 ). Support is also provided by the Smithso¬ 
nian Astrophysical Observatory, the MIT School of Sci¬ 
ence, the Raman Research Institute, the Australian Na¬ 
tional University, and the Victoria University of Welling¬ 
ton (via grant MED-E1799 from the New Zealand Min¬ 
istry of Economic Development and an IBM Shared Uni¬ 
versity Research Grant). The Australian Federal gov¬ 
ernment provides additional support via the Common¬ 
wealth Scientific and Industrial Research Organisation 
(CSIRO), National Collaborative Research Infrastruc¬ 
ture Strategy, Education Investment Fund, and the Aus¬ 
tralia India Strategic Research Fund, and Astronomy 
Australia Limited, under contract to Curtin University. 
We acknowledge the iVEC Petabyte Data Store, the Ini¬ 
tiative in Innovative Computing and the CUDA Center 
for Excellence sponsored by NVIDIA at Harvard Univer¬ 
sity, and the International Centre for Radio Astronomy 
Research (ICRAR), a Joint Venture of Curtin University 
and The University of Western Australia, funded by the 
Western Australian State government. 


APPENDIX 

A. DELAY TRANSFORM CONVENTIONS 



Local 

Meridian 


Figure Al. Radio interferometer delay conventions used in this 
paper. Two antennas labeled as Ai and A 2 are separated by vector 
b on the local tangent plane. The local meridian, local north and 
local east are shown for reference. Points labeled as ‘z’, ‘n’, ‘s’, ‘e’ 
and ‘w’ on the celestial sphere denote zenith, northward, south¬ 
ward, eastward and westward positions, respectively. r z , r n , r s , 
r e and r w denote the respective delays measured between Ai and 
A 2 . Throughout this paper, zenith is chosen as the phase center. 
Hence, r z = 0. If b is oriented eastward as shown, then r e > 0, 
r w < 0, and r n = r s = 0. Conversely, if b is oriented northward 
(not shown here), then r n > 0, r s < 0, and r e = r w = 0. 


B. COSMOLOGICAL H I POWER SPECTRUM 
Equation [T| can be equivalently expressed as: 

Vu(f) = jj A(a, /) I(a, /) Wi(f) e~ i2 ™° da, (Bl) 

sky 


where, s is measured with reference to a location on the 
sky referred to as the phase center , and u = (u,v,w) 
denotes the spatial frequency vector, w is aligned par¬ 
allel to the direction of the phase center, while u and v 
lie on the transverse plane perpendicular to it. For mea¬ 
surements that lie on this plane, we can choose w = 0 
without loss of generality and u effectively reduces to 
u = (u, u), a two-dimensional vector. Then, u is directly 
related to the transverse spatial wavenumber mode as: 


k ± = 


27 TU 

~cT ’ 


(B2) 


Figure IA1I illustrates the radio interferometer delays 
and conventions used in our paper, b is assumed to be 
on a coordinate system aligned with the local east, north 
(along local meridian) and upward (zenith) directions at 
the telescope site. Hence, a perfectly eastward oriented 
antenna spacing will observe objects in the eastern and 
the western skies at r > 0 and r < 0, respectively. Simi¬ 
larly, an object in the northern sky will appear at r > 0 
for an antenna spacing oriented northward. For all ob¬ 
servations used in this study, we use zenith as the phase 
center, for which r = 0 . 


where, D = D(z) is the transverse comoving distance at 
redshift z. 

Since we are concerned with a redshifted H I spectral 
line from cosmological distances, / is a measure of cos¬ 
mological distance along the line of sight. 77 , which is the 
Fourier transform dual of /, is used to denote the spatial 
frequency along the line of sight and has units of time. 
It is directly related to the line of sight wavenumber, 


2nr]j2iHoE(z) 

c( 1 + z) 2 


(B3) 














Foregrounds in EoR Power Spectra 


15 


where, /21 is the rest frame frequency of the 21 cm spin 
flip transition of H I, and Hq and E(z) = [Hm(1 + z) 3 + 
fik(l + z) 2 + Ha] 1 / 2 are standard terms in cosmology. 
This approximation holds under the assumption that the 
redshift range (or frequency band) is small enough within 
which cosmological evolution is negligible. Thus, 

Vu(v) = J Vu(f) W(f ) d/ (B4) 

represents the true spatial Fourier representation 
of the three-dimensional sky brightness distribu¬ 
tion. This approach h as been discussed in detail in 
iMorales fc Hewittl ([20041) . The spatial power spectrum 
of EoR H I distribution, P(fc_i_,fc ||), and V u {jj) a re re¬ 
lated by (IMorales fe Hewittl 120041 iMcQuinn et al.l [20061: 
iParsons et al.ll2012at) : 


P(k ± ,k {l )^\v u ( v )\ 2 


\X 2 AB J 


fD 2 AD\ 

\~AB~) 



(B5) 


where, A e is the effective area of the antenna, A B is the 
bandwidth, AD is the comoving depth along the line of 
sight corresponding to AF>, A is the wavelength of the 
band center, and k b is the Boltzmann constant. Thus, 
V u (rj) inferred from observations, in units of Jy Hz, can 
be converted into an equivalent cosmological H I power 
spectrum P(k |_, fey), in units of K 2 Mpc 3 or, more gener¬ 
ally, K 2 (Mpc/h) 3 , where h is the Hubble constant factor. 

Without loss of generality, the phase center can be as¬ 
sumed to be the zenith relative to the local tangent plane. 
Then u lies on this plane for measurements constrained 
to be on it. If the array of antennas are also coplanar 
lying on the local tangent plane, then u = 6/A. Under 
such circ umst anc es, e quations [T| and [2] closely resemble 
equations [Bl] and [B4] respectively. However, they are not 
quite identical to each other. It is because 6 is_indepen¬ 
dentoffreguency_while u is not. Par sons et all (|2012bl ) 
and iLiu et all (j2014al) have discussed the mathematical 
correspondence between the two. 
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